Statistical fluxes and the Curie- Weiss metal state 
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We predict a new state of matter in the triangular t-J model in a high doping regime. Due to 
the altered role of quantum statistics the spins are 'localized' in statistical Landau orbits, while the 
charge carriers form a Bose metal that feels the spins through random gauge fields. In contrast to 
the Fermi-liquid state, this state naturally exhibits a Curie- Weiss susceptibility, large thermopower, 
and linear-temperature resistivity, explaining the physics of Na^CoC^ at x > 0.5. A 'smoking gun' 
prediction for neutron scattering is presented. 



Introduction. — Among the great challenges in quan- 
tum matter physics is the question whether Mottness 
- the drastic change of Hilbert space structure due to 
strong local interactions - might form the condition for 
non Fermi liquid states of fermionic matter to occur— 
Among others, this is the point of departure of Ander- 
son's resonating valence bond (RVB) proposal for high- 
T c superconductivity 3 -. At the electron densities in the 
proximity of the half-filled Mott insulator the Mott pro- 
jections are most evident and it is much less evident 
whether they matter at all at very high(or dilute) elec- 
tron densities. Resting on a representational that makes 
explicit the altered nature of the quantum statistics due 
to Mottness we demonstrate here the existence of an in- 
ternally consistent, stable saddle point at high dopings 
on a triangular lattice. This is a spin-charge separated 
state where the doublons (carrying charge) and spinons 
(carrying spin) are coded as hard core bosons that com- 
municate with each other via statistical fluxes. The sur- 
prise is in the spin sector: the statistical gauge fields act 
like a uniform magnetic field causing Landau quantiza- 
tion of the spinon states. The spinons are localized in 
Landau orbits of the lowest Landau level (LLL) at low 
temperatures, with the effect that these behave as free 
Curie- Weiss spins with large entropy. The doublons get 
in turn scattered by the random gauge fluxes associated 
with the localized spinons, similarly as in gauge glass 
models. 

This might sound far fetched. However, the mi- 
croscopy matches that of the highly overdoped Na x CoC>2 
system^— with x > 0.5. A very strange metal is formed 
in this cobaltate, characterized by a high density of free 
Curie- Weiss spins that 'appear out of the blue' while it 
is a rather good hole-type thermoelectric material^ given 
its large thermopower combined with a relatively low re- 
sistivity. We will show here that the magnetic- (Fig. 1), 
thermoelectric- (Fig. 2) and resistive (Fig. 3) properties 
of the cobaltate are consistent with the state described in 
the previous paragraph. Moreover, we present a 'smok- 
ing gun' prediction that can be straightforwardly tested 
by experiment: the spin excitation spectrum should carry 
the fingerprints of the 'spontaneous' Landau quantization 
(Fig. 4). We notice that earlier attempts^— to explain 
the Curie- Weiss metal behavior of the cobaltates relied 
on the assumption of a very small bare bandwidth, while 



in our theory the small itineracy scale is emergent. 

Similar to the layered structure in the high-T c 
cuprates, the C0O2 layer is believed to play an essen- 
tial role in determining the low-energy physics in the 
Na^CoC^ compound. With electrons doping introduced 
by Na, the electron hopping at the partially filled t^g 
orbitals of the Co +3 ions and spin correlations between 
C0+ 4 ions in a C0O2 layer may be minimally described 
by the t-J modeU*~— on a triangular lattice with hopping 
integral t < 0l2r—. Here the Hilbert space is constrained 

by Yla c i<T C icr — 1j i- e -; eacn lattice site is either singly 
occupied by an electron (Co +4 ) or doubly occupied by 
electrons (Co +3 ) without allowing the empty site (Co +5 ). 
What we will be interested in is the highly overdoped 
regime of this model, where the RVBis— correlations 
induced by the superexchange coupling disappear and 
the hopping processes of both charge and spin become 
dominant. 

Sign structure. — To identify the new saddle point, 
we have to rely on a particular representation that is 
superficially reminiscent of the standard slave bosons 
but is actually quite different. This 'phase string' 
representation^ is making explicit the non Fermi-Dirac 
nature of the quantum statistics in doped Mott insula- 
tors. According to this formalism, the electron operator 
da can be fully 'bosonized' by c\ a = d\bi- a e~ l& i» E , 
where the doublon and spinon creation operators, d\ 



and b\ a . are both bosonic, which satisfy an equality 

nj ff = 1 with nf and n\ a 



(Mott) constraint nf + £ CT n" t . 



denoting the doublon and spinon occupation numbers, 
respectively. Here the fermionic commutations of the 
electron operators are ensured by the topological phase 
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(zi is the complex coordinate of site i). In this represen- 
tation, the hopping and superexchange terms of the t-J 
model become 
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In this way, the remnant sign structure after Mott pro- 
jection has been made explicit while it is precisely rep- 
resented by the topological link variables, Afj and Af- 
- Obviously there would be no 'sign problem' should 
Afj = Afj = in such a fully bosonized model with 
t < 0. Here A" 4 and Af .■ satisfy 



ij 
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E 4 = -E(W) 
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(3) 
(4) 



where dS* denotes the boundary of an area S. 

Electron fractionalization saddle-point. — The precise 
sign structure identified above will be crucial in con- 
structing the following saddle-point which respects the 
gauge invariance associated with Afj and Afj, together 
with time-reversal and spin rotational symmetries. Since 
the RVB pairing is irrelevant at high doping, the gauge- 
invariant Dij and Bfj in Eqs. ([1]) and © will be nat- 
ural order parameters, resulting in an effective H c g = 
Ha + H s : 



exp. 
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H d = -t d E e " <0y (e iA v4dj)+H.c. (5) 
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H s = -t, e ~ iaii (e- ia <blb ja )+H.c. (6) 
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FIG. 1: The spin part of the present electron fractionaliza- 
tion state exhibits a Curie T-dependent uniform susceptibil- 
ity Xu = C/T (full squares) at x = 0.69. The Weiss term 
in x« appears as an RPA correction (full circles). The Fermi 
liquid (dotted) and 'classical limit' (dashed) results are also 
shown for comparison. Inset: the coefficient C (solid) and 
the 'classical limit' C c \ (dashed) vs. doping x, in which the 
experimental data are marked by full triangles^ ' 21 ' 22 . 



where t d = —tB > 0, t s = —tD + JB/4 > 0. Here oy chemical potential). As clearly illustrated in Fig. Q]by 

with 6 = 



represents the U(l) gauge fluctuations around the saddle- 



point: ~ De mi i and Bfj ~ (B/2)e to «. It can be 
showr*i£ that the transverse fluctuations of ay will get 
suppressed^ to ensure the stability of the saddle-point 
and the flux binding conditions Eqs. ([3]) and ((4]), imply- 
ing spin-charge separation. Note that the fractionalizcd 
state as governed by H e g at the gauge fixing ay = is 
still an intrinsic gauge model involving the mutual Chern- 
Simons gauge fields, Afj and Afj. 

Mean-field approximation. — The spin dynamics is gov- 
erned by H s , where Afj given in Eq. (@| is only density- 
dependent and can therefore to leading order approx- 
imation be treated as a smearing flux of 7r(l — x) per 
unit cell. The resulting 'mean-field' state of H s is Lan- 
dau quantized with spins being statistically 'localized' in 
the cyclotron orbits, in sharp contrast with a degener- 
ate Fermi liquid state or a fully localized classical state 
at t — > 0. In particular, the characteristic bandwidth 
vanishes, while t remains relatively large, when all the 
spinons stay in the degenerate LLLs at low temperatures, 
rendering a systematic 'scaling' behavior as will be shown 
bclow^S. 

Curie-Weiss uniform susceptibility. — One peculiar 
property is exhibited by the uniform spin susceptibil- 
ity, Xu = ^7r^I]m( n ™ + l ) n m. where = k B T and 
n m = l/[e^( £m ~^ s ) — 1] is the Bose distribution for 
spinons at state m with energy E m obtained by the afore- 
mentioned mean-field solution based on H s (fj, s is the 



1/Xui h follows a Curie- Weiss law Xu - 
(full squares). A finite Weiss term O ~ 3(1 — x)J/k-Q 
(full circles) is generated by including an RPA correction 
from Hj (by fitting with the experimental 0£, we esti- 
mate J ~ 70 K). The coefficient C is x-dependent, and 
at T — > one finds C = 2(1 — x)^/Hb- As shown in the 
inset of Fig. [I] C is in quantitative agreement with the 
experimental data £ 21 ' 22 in the Curie- Weiss regime of the 
cobaltates, which is independent of other parameters in 
the model, like t and J. 

By contrast, a Pauli-like susceptibility is expected for 
the Fermi-liquid state (dotted line in Fig. [TJ obtained 
with a bare t = — O.leV— ). It is particularly instructive 
to compare this with the classical limit of the t- J model at 
t (J = 0), where all the electrons are fully localized 
as (1 — x) N free moments, contributing to a Curie's law 
X^ 1 = C c i/T (dashed line in Fig. [TJ. But one finds C d = 
C/2, i.e., only about the half of the values of both the 
experimental and the present theory (cf. the inset of Fig. 
[1]). Clearly the peculiar quantum effect of the present 
bosonic spinons is responsible for the enhancement of C 
from C c i at low T (only at high-T limit, the 'classical' 
Xu can be recovered as n m — > in Xu)- 

Thermopower. — In Fig. [5] the thermopowers predicted 
by the Fermi-liquid state as well as the 'classical limit' 
of the t-J model are shown by dotted and dashed curves, 
respectively In the latter case, with the bandwidth van- 
ishing, the thermopower reduces to the so-called Hcikcs 
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formula: Q c \ = (fce/e) In rp^l wn i cn i s proportional to 
the entropy per electro n 25 ! 26 . Both deviate strongly from 
the experimental result (full triangles in Fig(2J) in oppo- 
site ways. 

In the present saddle-point state, the thermopower will 
satisfy the Ioffc-Larkin combination rule Q = Qd + Qs ac- 
cording to Eqs. and ©, with the spinon contribution 
Q s dominant over the doublon part Qd- Here a 'large' 
Q s originates in the degeneracy of the LLLs. This leads 
to a Heikes-like formula Q s = —fx s /eT, associated with 
the spinon entropy, which follows given that the energy- 
particle current correlator S 12 = at low temperature in 
the Kubo formula Q 3 = -^frr - Note that the 

above Heikes-like formula is still valid even when S 12 be- 
comes finite, for instance when the spinons are excited to 
the higher Landau levels at high-T, because the current- 
current correlator S 11 — > oo 18 ' 19 . 

The calculated Q s is presented in Fig. [2] (solid- 
squares). This saturates at low-T to a 'universal' value 
Q s0 = In 2 p» 60 pN/K due to the artifact that the 
bandwidth of the LLLs remains zero. But a weak fluc- 
tuation of the gauge flux of Afj around the mean- value 
can easily cause a broadening lifting the exact LLL de- 
generacy, resulting in a vanishing Q s at T — > 0. A typical 
example is indicated by the full circles in Fig. [2j with Ay 
simulated by a random flux of [— 0.05-7T, 0.057r] around the 
mean flux of 7r(l — x) per unit cell. One finds a quanti- 
tative agreement with the experimental measurement by 
taking t s = 100 K. Note that Qd contributed by dou- 
blons is not included here as it is usually much weaker 
in magnitude with a quite flat T-dependent as calculated 
by the Kubo formula^. 

To understand more clearly the origin of Q s from the 
spin degrees of freedom, one may apply a strong in-plane 
magnetic field H to polarize (freeze) the spins via Zeeman 
energy, as already accomplished experimentally^. Denote 
Q s (oo) as the contribution from the remaining configu- 
rational entropy at H — > oo and define 

Q S (H)-Q _ 

Q s (0)-Q s (oo) ■ {! 

Then a is found to be a well-defined scaling curve of 
H/T, as shown by the solid curve in the inset of Fig. 
[21 which fits well with the pure spin entropy of free 
moments: [In [2 cosh(w)] — u tanh(it)] / In 2, where u = 
(2.2/j.bH) I (2ksT) in excellent agreement with the ex- 
perimental data as given in Ref£. 

Resistivity. — The resistivity satisfies the Ioffc-Larkin 
combination rule p = pd + p s in which p s = due to 
the aforementioned "Meissner effect" respons o 18 ' 19 in the 
spinon subsystem such that the doublon subsystem will 
solely contribute to p. Governed by Hd, the doublons will 
experience the novel scattering generated from the gauge 
potential Afj, which describes ±7r flux-tubes bound to 
the spinons with a = ±. In the low-T regime where 
the spinons are in the degenerate Landau orbits, the cor- 
responding flux-tubes seen by the holons will distribute 




FIG. 2: Thermopower Q s contributed by the spin part at 
x = 0.71. Full squares: the mean-field solution; Full circles: 
the effect of doublon density fluctuations is considered (see 
text). The experimental result (full triangles) at the same 
x as well as the Fermi liquid (dotted) and 'classical limit' 
(dashed) results are shown for comparison. Inset: The scaling 
curve (solid) for a defined in Eq. 10 vs. H/T and the spin- 
entropy of free moments (dashed) , which well accounts for the 
experiment measurement. 




T(t d ) 

FIG. 3: The resistivity is contributed by doublons in the 
presence of random flux tubes ±7r attached to spinons in the 
LLLs, which numerically treated as a random distribution of 
±(1 — S)tv fluxes per unit cell. 



randomly in space. Based on a numerical calculation 
where A|- is treated as randomly distributed ±(1 — x)ir 
fluxes per unit cell, we obtain a metallic behavior of 
Pd with a lincar-T dependence over a large temperature 
regime as shown in Fig. G2 originating in a scattering rate 
similar to a case previously studied 2 ^ in a high-T regime 
in the context of high-T c cuprates. 

Predictions and discussion. — The fractionalized 
saddle-point state governed by (O and © has been 
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FIG. 4: Dynamic spin structure factor S(q, uj) predicted by 
the mean-field theory at x — 0.75. Right panel: the Lan- 
dau level structure exhibited in the local (q-integrated) dy- 
namic structure factor Sl(w); Left panel: the corresponding 
q-dependence of S^q, cu) for the first three energy levels. 



shown to gcnerically exhibit a systematic scaling behav- 
ior: p oc T, Xu oc 1/T, a large thermopower Q s ~ — In 2 
as well as the scaling law in a = a(H/T), associated 
with the peculiar Landau quantization effect in H s . Our 
theory allows us to make one further strong prediction: 
the spin sector Landau levels can be in fact directly 
observed by inelastic neutron scattering! The dynamic 
spin structure factor 5(q, w) can be easily computed 
assuming the straightforward Landau quantization and 
it should look like Fig. [4j a tower of rather narrow iron 
dispersive bands of spin fluctuation (right panel) , with a 
momentum dependence as indicated in the left panel. 

Another interesting consequence of the spin-charge 
separation is in the figure of merit, used to quantify the 
efficiency of a thermoelectric device. This can be written 
in the following form 

ZT ee 9^- = (Qs + Qd ^l Ld (8) 
pn 1 + K p h/nd 



since the spinons do not contribute to the thermal con- 
ductivity and resistivity (k s = and p s = 0)i&. Then a 
large spinon thermopower Q s (>> \Qd\) will play an im- 
portant role, independent of the Lorenz number of the 
doublons La = T as well as the doublon (ka) and 

phonon (K p h) thermal conductivities. We findi^ Ld to be 
one order of magnitude smaller than the Lorenz number 
of a Fermi liquid. This is in fact consistent with exact 
diagonalization result on the t-J models. The combina- 
tion of these factors may therefore lead potentially to an 
exceptionally large Z. 

Finally, we already alluded to the fact that various 
factors leading to a broadening of the Landau levels are 
ignored in the present mean-field treatment^. Accord- 
ingly, the physics we have discussed here is associated 
with an intermediate temperature regime, bounded by 
the Landau gap scale at high temperature and the onset 
of further cooperative phenomena at low temperatures. 
In the latter regime the specific heat measurements have 
also show n 29 ' 30 an anomalous increase, presumably due 
to the afore-discussed entropy enhancement which, how- 
ever, will be sensitive to the detailed Landau level broad- 
ening. Eventually the spin degeneracies associated with 
the Landau states will be lifted, and this might be the 
origin of some new emergent orders at very low temper- 
atures. One infers from the momentum dependence of 
5(q, oj) at the LLLs (Fig. 2]) that the in-plane ferromag- 
netic correlation length will increase with x, such that a 
small interlayer superexchange J^may eventually drive 
an in-plane ferromagnetic order with antiferromagnetic 
ordering along different layers, explaining the A- type an- 
tiferromagnetic order observed by neutron scattering at 
x > 3/421 The bosonic doublons may also condense 
at a sufficiently low-T to make the system more Fermi- 
liquid like, reconciling with the Wiedemann-Franz law 
observed in Nao.7Co02 at T < 1 However, our un- 
conventional 'building material' also leaves room for less 
conventional forms of order, including uncommon charge 
orders that will be discussed in detail elsewhere. 
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